Interaction-Induced Localization of an Impurity in a Trapped Bose Condensate 
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We study the ground state properties of a trapped Bose condensate with a neutral impurity. By 
varying the strength of the attractive atom-impurity interactions the degree of localization of the 
impurity at the trap center can be controlled. As the impurity becomes more strongly localized the 
peak condensate density, which can be monitored experimentally, grows markedly. For strong enough 
attraction, the impurity can make the condensate unstable by strongly deforming the atom density 
in the neighborhood of the impurity. This "collapse" can possibly be investigated in bosenova-type 
experiments. 
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I. INTRODUCTION 

The study of impurities immersed in liquids and solids 
has a long history. In 1933, Landau predicted, using 
quantum mechanical arguments, that the localization of 
electron impurities in a crystal could be used to probe the 
activation energy of solids p) . Electron impurities have 
also played a key role in the study of liquids, in particular 
liquid 4 He Q. More recently, the study of doped meso- 
scopic helium clusters has attracted much attention 
Some atom impurities reside on the cluster surface while 
others migrate to the center of the helium cluster. Spec- 
troscopic measurements of molecules located at the cen- 
ter of the cluster have, e.g., shown unambiguously that 
4 He clusters with about 60 atoms are superfluid |g. 

Recently, the study of impurities immersed in a 
gaseous, coherent atom background has become possi- 
ble 0, @. Theoretical studies on, e.g., ion impurities 
in a condensate have been initiated |a |9( , raising ques- 
tions about the appropriate treatment of systems with 
long-range interactions (unlike short-ranged atom-atom 
potentials, which behave as 1/r 6 for large interparticle 
distances, atom-ion potentials fall off as 1/r 4 ). Here, we 
consider a neutral impurity in an inhomogeneous Bose 
gas, assuming contact atom- impurity interactions |10| . 
Treatments for more complicated atom-impurity inter- 
actions exist 0, 0, 0; the results may, however, be 
model-dependent. Our self-consistent mean-field treat- 
ment provides a first step towards a systematic under- 
standing of impurities in a Bose condensate. We also 
discuss a simple variational treatment which reproduces 
the key features of the self-consistent results. We point 
towards possible experimental signatures of our predic- 
tions, which will be aided by the possibility of tuning the 
atom-atom and atom-impurity interactions in the vicin- 
ity of a Feshbach resonance by application of an exter- 
nal magnetic field [13, H^. This tunability is unique to 
gaseous condensate-impurity systems; it does not, for ex- 
ample, exist in helium where the interaction strength is 
set by nature. 

We consider a weakly-interacting Bose condensate in 
a harmonic trap, doped with a single impurity. For now, 
we assume that the impurity feels no external trapping 
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FIG. 1: Phase diagram for trapped Bose gas with a sin- 
gle impurity, which feels no confining potential, as a function 
of the number of atoms N and the atom-impurity scattering 
length a a i for equal atom and impurity mass, i.e., mi = m a . 
The phase diagram contains three regions: in region (A) the 
impurity is unbound; in region (B) the impurity is localized 
(the localization is "weak" for comparitively small \a a i\ and 
"strong" for comparatively large \a a i\, see Sec.|n]for details); 
and in region (C) short-ranged physics becomes relevant. Re- 
gions (A) and (B) are separated by a critical value a i, c l (up- 
per bold solid line), which is approximately independent of 
a aa - Regions (B) and (C) are separated by a critical value 
a a i,c2, which is shown for a aa — 0.005a ho (lower bold solid 
line), a aa = (dotted line), and a aa ~ 0.05ah o (dashed line). 



potential; later, we discuss how the presence of an im- 
purity trapping potential modifies the results. Figure 1 
shows the equilibrium "phase diagram" |l6j | determined 
within mean-field theory as a function of the number 
of atoms N and the atom-impurity scattering length 
a a i- The phase diagram separates into three distinct re- 
gions: (A) For a a i > a ai cl , the impurity is unbound and 
can move away from the trapped atom cloud. (B) For 
Oai.ci > o-ai > dai.r.2-, the impurity is localized, i.e., bound 
to the atom cloud |T3 - (C) For a ai < a a , Lc2 , short-range 
physics, which cannot be described within mean-field the- 
ory, becomes relevant. Regions (A) and (B) are separated 
by a TV-dependent critical value a a i, c \ (upper solid bold 
line in Fig. 1), which is approximately independent of 
the atom-atom scattering length a aa . Since the impurity 
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feels no trapping potential, interaction-induced localiza- 
tion of the impurity occurs only if a ai is more attrac- 
tive than a a i,d- Regions (B) and (C) are separated by 
a iV-dependent critical value a Q i jC 2, which also depends 
on the atom-atom scattering length a aa . The lower bold 
solid line in Fig. 1 shows a a i, C 2 for a aa — 0.005a/ lo , the 
dotted line that for a aa = 0, and the dashed line that 
for a aa = 0.05aho- If o-ai is more negative than a a i, C 2, 
the attractive atom-impurity interactions can "collapse" 
the condensate, pulling atoms into a short-ranged state 
about the impurity. 

The next section outlines the self-consistent mean-field 
treatment used to calculate the phase diagram shown in 
Fig. 1. Section ITTT1 develops a simple variational frame- 
work, which reproduces the key features of the full self- 
consistent mean-field treatment. Finally, Sec. II VI dis- 
cusses possible experimental realizations of the systems 
under study and concludes. 



II. SELF-CONSISTENT MEAN-FIELD 
TREATMENT 

We describe N atoms of mass m a in the presence of a 
harmonic trapping potential with angular frequency LOho 
and a single impurity of mass m,, which feels no exter- 
nal potential, within mean-field theory. Assuming that 
the atom-atom and atom-impurity interactions can be 
described by contact potentials, the many-body Hamil- 
tonian reads 
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where U qp — 2ith 2 a qp /m qp , m qp = m q m p /{m q +m p ), and 
(q,p) — (a, a) or (a,i). In Eq. QJ, Xj and x*i denote the 
position vectors of the jth atom and the impurity, respec- 
tively. We approximate the ground state wave function 
$ as a product of single-particle wavefunctions, 
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and derive a set of coupled Hartree-Fock equations, 
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FIG. 2: (color online) Self-consistent mean-field results ob- 
tained for attractive atom-impurity interactions, iV = 10 4 , 
a aa = 0.005a,ho, and m, = m a . Panel (a) shows the chem- 
ical potentials e; and e a as a function of a a i- Solid lines 
in panels (b) and (c) show the effective potentials V e ff,i(r) 
and V e ff,a(r), respectively, for a few selected atom- impurity 
scattering lengths (see legend). Dashed lines show the corre- 
sponding chemical potentials [note that the change of e a is not 
visible on the scale chosen in panel (c)]. The insets of panels 
(b) and (c) show the corresponding impurity wave function 
ipi(r) and atom density |^ a (r)| 2 , respectively. The critical 



values a a i, c i and a a i, C 2 of this system are 
and — 0.062ah o , respectively. 
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Here, aho denotes the oscillator length (a/ lo = yJh/m a ujho 
and r = ah x), and e a and the chemical potentials 
(or "orbital energies") of the atoms and the impurity. 
The coupled mean-field equations are equivalent to those 
for a two-component condensate |l8j if one replaces one 
of the two components by a single impurity. In writ- 
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ing Eqs. J3J and (@J, we have implied spherical symme- 
try, ^ a .i{r) = VwO) = with Jo°° \Xa,i{r)\ 2 dr = 1. 
For repulsive atom-impurity interactions, not considered 
here, symmetry-breaking states can exist |Tflj . 

The impurity feels an effective potential V e ff % [defined 
as the second term in square brackets on the left hand 
side (LHS) of Eq. QJ], which is created by the atom 
density \ip a \ 2 - The impurity density \ipi\ 2 enters Eq. (J3J 
and creates, together with the trapping potential and 
the atom density itself, an effective potential V4//, a [de- 
fined as the last three terms in square brackets on the 
LHS of Eq. J2J]. For weak atom-impurity interactions, 
the condensate atoms act to a good approximation as 
a static background with which the impurity interacts. 
However, as the strength of the atom-impurity interac- 
tions increases the full coupled nature of Eqs. © and Q 
becomes important. 

We discuss the self-consistent solutions to Eqs. (0 and 
Q, obtained by numerical means, for a specific set of 
parameters. The behavior is qualititatively similar for 
other parameters. Figure |2Ja) shows the chemical po- 
tentials €i and e a as a function of a a i for N = 10 4 and 
&aa = 0.005a^ o , and m s = m a . Equal atom and impurity 
masses can be realized experimentally by, e.g., promoting 
a single condensate atom to a different hyperfine state [6j . 
Figure |2Ja) shows that the chemical potentials Cj and e a 
change approximately linearly with the atom-impurity 
scattering length for —0.005 > a ai /aho > —0.05. This 
linear behavior is what one would expect from a per- 
turbative treatment. To visualize the system's behavior, 
Figs.EJb) arid (c) show the effective potentials 
and V e f f t a (t) for a few selected atom- impurity scattering 
lengths. Figure 2(b) shows that V e ff,i becomes deeper 
as dai/aho goes from —0.01 to —0.03 to —0.05. Accord- 
ingly, the impurity wave functions ipi, shown in the in- 
set of Fig. 2(b), become more localized as \a a i \ increases. 
Although \tpi\ 2 changes significantly as a a i/aho goes from 
—0.01 to —0.05, V e ff a and IV'al 2 change only slightly [see 
Fig. He)]. 

For a a i/a,ho —0.05, the impurity chemical poten- 
tial ei (and, to a lesser degree, the atom chemical poten- 
tial e a ) changes in a non- linear, i.e., non-perturbative, 
fashion. For the parameters at play here, this de- 
fines the regime of strong atom-impurity coupling. To 
highlight the dramatic changes of the system in this 
strongly-coupled regime, Figs. 2(b) and (c) show self- 
consistent effective potentials for three nearly identi- 
cal atom-impurity scattering lengths, i.e., a a i/aho = 
-0.06,-0.061,-0.0615. The peak impurity density 
grows with increasing \a a i\ and creates a "hole" at the 
center of V e f , which in turn causes the atom density to 
grow a "bump" at the trap center [see inset of Fig. 2(c)] 
with a length scale of roughly the condensate healing 
length £. The healing length is given by the competition 
between the kinetic energy and the condensate's mean- 
field energy, £ = l/\/8np a a aa [2?i|. where p a denotes the 
peak density of the atoms, p a — N\t/j a (r — 0)| 2 . For 
the parameters of Fig. [21 £ ~ 0.26a,h o - Since the healing 



length £ is the scale over which the condensate "reacts" 
to spatial perturbations, it is natural that the atom den- 
sity develops a variation near the trap center of size £ . 

The inset of Fig.EIc) illustrates the peak atom density 
growth with increasing |a a j|. To quantify this growth, 
we calculate the excess number of atoms AN associated 
with the bump of the atom density. In analogy to a 
homogeneous system we define AN as 

AA = AttN f °[|<^°(r)| 2 - h£" =0 (r)l V*", (5) 
Jo 

where the atom wave function 4>a ai ^°{ r ) 1S calculated 
self-consistcntly for a Bose gas with finite atom-impurity 
scattering length and ipa ai= °( r ) f° r a Bose gas with van- 
ishing atom-impurity scattering length (for the same A 
and m a ). When evaluating Eq. (J5J for a specific system, 
we choose the cutoff radius r c to roughly coincide with 
the r-value at which the bump of the atom density starts 
growing. Triangles in Fig.EJshow the resulting number of 
excess atoms AN for N — 10 4 , m t = m a , a aa — 0.005a/j O 
(the same parameters as in Fig- EJ and r c — la/ lo for dif- 
ferent values of the atom- impurity scattering length a a %- 
The number of excess atoms increases roughly linearly 
with increasing |<z j|. Just before the onset of instabil- 
ity at a j )C 2 ~ — 0.062a^ o , the number of excess atoms 
AA reaches 12, which corresponds to 0.12 % of the to- 
tal number of atoms. For comparison, a dashed line in 
Fig- El shows an estimate for the number of excess atoms 
derived for a weakly-interacting impurity-doped homoge- 
neous Bose gas 0, 

AN = - — — . (6) 

Figure El shows good agreement between the number 
of excess atoms A7V calculated for the inhomogeneous 
impurity-doped condensate (triangles) and the analytical 
expression, Eq. |JBJ|. This suggests that Eq. © describes 
the number of excess atoms for large enough, weakly- 
interacting inhomogeneous condensates quite accurately. 

Finally, if a a t becomes more negative than a critical 
value of a a i,c2 ~ — 0.062a^ o , we no longer find a self- 
consistent solution to Eqs. © and Q. This implies that 
the condensate collapses, i.e., atoms are drawn into a 
short-ranged state about the impurity. It appears likely 
that this collapse involves only a fraction of the conden- 
sate atoms, but a definite answer lies beyond the scope 
of the present work. Just as in the case of pure atomic 
condensates with negative atom-atom scattering length, 
mean-field theory predicts the onset of collapse for our 
coupled equations but cannot describe the system's be- 
havior in the regime dominated by short-range physics. 

We now estimate the critical value a a i,ci) which sep- 
arates the unbound impurity phase from the localized 
impurity phase for n%i = m a (see upper bold solid line 
in Fig. 1). The impurity feels a strictly short-ranged 
potential, i.e., V e ff,i[r) falls off faster than a power- 
law of r. Since the impurity equation, Eq. Q, is lin- 
ear, we can compare the volume-integrated strength of 
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calculated, Eq. (5) 
— estimate, Eq. (6) 
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FIG. 3: Triangles show the number of excess atoms AiV 
calculated from Eq. © using the self-consistent mean-field 
solutions for N = 10 4 , a aa = 0.005ah o and m a — rrii (the 
same parameters as in Fig. 2), as a function of the atom- 
impurity scattering length a ai . For comparison, a dashed line 
shows the analytical estimate, Eq. JSJ. 



V e ff t i with the corresponding critical value, 
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for forming a bound state in three dimensions in a po- 
tential of range b. If we identify the range b with a^ Q , 
we find that the critical value of a a i scales as 1/N, i.e., 

a a i,d(N)/a ho = ~2iiv ~ ~Hr- For Positive a aa , the 
atom cloud is somewhat larger, and our estimate will be 
off by a numerical factor of order 1. The upper bold 
solid line in Fig. shows our analytical estimate. The 
results of our numerical calculations are consistent with 
this analytical estimate. The critical value a a j )C i might 
be difficult to observe experimentally since the transition 
from region (A) to region (B) involves a diverging length 
scale. Furthermore, it might be difficult to experimen- 
tally realize a trapping setup with tunable atom-impurity 
scattering length a a i for which the impurity feels no con- 
fining potential (see also Sec. HVjl . 

The behavior of the impurity-doped condensate was 
illustrated in Fig. 2 for N = 10 4 atoms. We find sim- 
ilar qualitative features for a smaller number of atoms, 
including the disappearance of the mean- field solutions. 
The critical values a a i, c i and a Q i,c2 vary with N as shown 
in Fig. 1. We note that the notion of a condensate heal- 
ing length, and thus the discussion of the disappearance 
of the mean-field solution at the point when the impurity 
becomes more tightly localized than this scale, becomes 
less meaningful for small enough number of atoms. It 
will be interesting to further investigate the properties of 
an impurity immersed in a small condensate since such 
systems can be realized experimentally with the aid of 
optical lattices. 



III. VARIATIONAL TREATMENT 

To further illustrate how the neutral impurity alters 
the atom cloud, we minimize the total energy variation- 
ally for the Hamiltonian given in Eq. Q and the wave 



function given in Eq. © with 

ipi cx e~ pir \ and ip a oc {e~ p ^ +ce" ftr! ). (7) 

The variational parameters pi and p a determine respec- 
tively the width of the impurity and of the atom wave 
function. To be able to describe the growth of the atom 
peak density in the strongly-interacting regime and the 
collapse of the condensate, ipa contains an additional 
gaussian with two more variational parameters, the rela- 
tive amplitude c and parameter pj,. We restrict pb to be 
greater than p a to separate the background condensate 
cloud from the more localized condensate bump. 

Figure 0] shows the results of our variational calcula- 
tions for N = 10 4 , a aa — 0, and to$ = m a as a function of 
a a i- To reduce the parameter space we set p a = 0.5aT o ; 
we checked that allowing p a to vary changes its value 
only little. The optimal values of the remaining three 
variational parameters pi, pb and c are shown in pan- 
els (c) and (d) of Fig. 4 by triangles, and the corre- 
sponding chemical potentials of the impurity and atom, 
respectively, are shown in panels (a) and (b) by trian- 
gles. The variational analysis predicts a critical value 
of a<M,c2 ~ — 0.020a/j O ; at this critical value a a i jC 2, the 
local minimum in the variational energy disappears as 
the variational energy becomes unbounded from below. 
For comparison, our self-consistent solutions to Eqs. © 
and which are shown in panels (a) and (b) by cir- 
cles for comparison, predict a somewhat less attractive 
critical value, i.e., a a i^2 ~ — 0.016a/j O . This is to be ex- 
pected since the self-consistent total energy (not shown) 
provides a better lower bound than the variational en- 
ergy. The variational parameters pi, pb and c shown in 
panels (c) and (d) of Fig. 4 nicely illustrate the degree 
of impurity localization. For small |a a j|, the amplitude c 
is negligible, indicating that the presence of the impurity 
barely affects the condensate. As \a a i\ increases, the im- 
purity becomes more tightly localized, i.e., pi increases 
(note that the width of the impurity density scales as 
1/ s/Pi), which in turn drives the growth and localization 
of the condensate bump, i.e., c andpb also increase. Since 
Pi drives the increase of pb, Pb necessarily increases slower 
than pi with increasing |<z a j|. 

To connect our results for the variational parameter 
c with the full self-consistent solutions, we fit our solu- 
tions to Eqs. and (0} to the wavefunctions of Eq. (J7J 
with the proper normalization, treating p.;, pb, and c as 
fitting parameters. The circles plotted in Fig. 4(d) show 
the resulting values of c extracted from the self-consistent 
solution. To a very good approximation, c 2 describes the 
percentage change in the peak condensate density for the 
system with non- vanishing a a i as compared to the system 
with vanishing a a i (assuming we keep a aa and the number 
of atoms N fixed). Figure 4(d) shows that, just before 
collapse at a ai:C2 , c w 0.53 for the full self-consistent so- 
lution to Eqs. © and and c w 0.68 for the variational 
solution. These values of c correspond to changes in the 
peak condensate density, as compared to the condensate 
without impurity, of greater than 25%. We note that a 
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FIG. 4: (color online) Triangles show the chemical potentials 
ei and e a [panels (a) and (b)], and the parameters pi,, pi and 
c [panels (c) and (d)] for N = 10 4 , a aa = and rrii = m a 
as a function of a a i, obtained from the variational treatment. 
For comparison, circles in panels (a) and (b) show £j and e a , 
and those in panel (d) the parameter c obtained by fitting the 
self-consistent solutions of Eqs. and @ ( see text). 



similar growth of the peak condensate density is seen in 
the inset of Fig. 2(c) for the same number of atoms, i.e., 
N = 10 4 , but non-vanishing atom-atom interactions, i.e., 
a aa = 0.05a ho . 

The variational wavefunction given in Eq. J7J is best- 
suited to describe the case of a aa = 0. For a aa ^0 
the atom cloud deviates from a gaussian, and for strong 
enough interactions a gaussian form for the atom cloud is 
a poor approximation. Consequently, as the atom-atom 
interactions increase, we find that the simplistic varia- 
tional wave function given in Eq. J7J cannot describe the 
tightly-localized impurity at the trap center prior to the 
onset of collapse. Nonetheless, for the case a aa = dis- 
cussed above, the variational treatment reproduces the 
key features of the full self-consistent solution and pro- 
vides us with further insights. In particular, the form of 
the variational wavefunction, Eq. Q, is useful in visual- 
izing how the condensate develops features characterized 
by a length scale much smaller than the oscillator length. 
Furthermore, the disappearance of the local minimum as 



the variational energy becomes unbounded from below is 
another indication, along with the disappearance of the 
self-consistent solutions, of the collapse of the conden- 
sate. 



IV. DISCUSSION AND CONCLUSION 

SectionslrTland llTll discuss the behaviors of a single neu- 
tral impurity, which feels no external confining potential, 
immersed in a trapped condensate. If the impurity feels 
an external trapping potential with angular frequency 
u>i, which might be the case in an experiment (see be- 
low), region (A) in Fig. 1 is absent, i.e., the impurity is 
always localized due to the presence of the external po- 
tential. Assuming that the impurity trapping potential 
has a characteristic length which is larger than roughly 
the condensate healing length (for large enough number 
of atoms), the comparatively strong impurity localiza- 
tion prior to collapse and the crossover from region (B) 
to region (C) in Fig. 1 are, however, nearly unaltered. 
For example, for N — 10 4 , a aa = 0.005a^ o , rrii = m a 
and us a = U3%, the onset of collapse occurs at the same 
critical value of a a i lC 2 ~ —0.0620^0 that we find without 
impurity trap. 

Impurity-doped condensates can be realized experi- 
mentally with present-day technology [||. If one con- 
siders a magnetically trapped condensate, an impurity 
can, e.g., be created by promoting one of the condensate 
atoms to a different hyperfine state. The promoted atom 
may or may not feel the magnetic confinement. Alter- 
natively, one could implant a different atom, magnetic 
or non-magnetic, into the cloud. Such systems have the 
disadvantage that the atom-impurity interactions can- 
not be tuned via a magnetic Feshbach resonance. To 
take advantage of the tunability of interspecies scatter- 
ing lengths |22j, one can consider an optical potential 
red-detuned with respect to the atoms and the impurity. 
In such an experimental realization both the atom and 
the impurity would feel trapping potentials. 

As the atom-impurity interactions are tuned closer to 
Qai,c2, the growth of the peak atom density at the center 
of the trap can potentially be monitored experimentally 
in expansion experiments. Since the condensate bump 
at the trap center involves only a few atoms, direct de- 
tection of the changes in the peak density may, however, 
be non-trivial. We suggest that the impurity-doped con- 
densate could alternatively be probed in a bosenova-type 
experiment which applies a sequence of time-dependent 
magnetic field ramps p^ . By tuning the atom- impurity 
scattering length to a large negative value, one could ex- 
perimentally induce collapse and consequently density os- 
cillations, which might involve a significant fraction of the 
condensate atoms. 

A key result of our study is that the degree of localiza- 
tion of the impurity at the trap center in region (B) of 
the phase diagram (see Figs. ^ and [5J can be controlled 
by varying the atom-impurity scattering length, i.e., the 
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width of the impurity wave function for a system with 
m a = rrii can be varied from a size much greater than the 
oscillator length a^ to a size significantly smaller than 
dho- I n addition to changing the atom- impurity scatter- 
ing length, one can consider unequal atom and impurity 
masses, e.g., a Cs atom immersed in a Na condensate. 
Not surprisingly, as the impurity mass increases, the de- 
gree of impurity localization also increases. The localized 
impurity itself may present the possibility for forming 
interesting single-atom devices, perhaps using the impu- 
rity's spin degrees of freedom. The favorable coherence 
properties of Bose condensates may make the localized 
impurity states viable for quantum computing schemes. 
Furthermore, extensions to two or more impurities will 
allow one to consider the role of condensate-mediated in- 



teractions between impurities. 

Finally, we return to our finding that a single neu- 
tral impurity can deform the condensate sufficiently to 
induce a collapse which may only involve a fraction of 
the condensate atoms. The resulting collapsed state may 
be related to the mesoscopic droplets that have been pre- 
dicted to form about an ion immersed in a condensate • 
More work is needed to fully understand these ion states, 
the collapsed states predicted in the present work, and 
possible connections between the two. 

We gratefully acknowledge fruitful discussions with A. 
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mans, hospitality of the Institute for Nuclear Theory at 
the University of Washington, and support by the NSF, 
grant PHY-0331529. 
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